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symmetry breaking parameter (magnetic flux). The applicability of the obtained 
distributions for the spectral statistic of disordered metals is discussed. The most 
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general N x N random matrix model. Furthermore, we find an unusual logarithmic 
behavior of the two point current-current correlation function at small flux. 
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I. INTRODUCTION 



The statistical properties of a single electron energy spectra in disordered mesoscopic 
systems have been the subject of intensive study during the past decade. Apart from a 
fundamental interest on it's own, the spectral statistic is closely related to such observable 
characteristics as persistent currents, anomalous magnetization, and conductance fluctua- 
tions. After the seminal studies of Efetov ||T| and Altshuler and Shklovskii it became clear 
that the spectral statistic in the diffusive regime, where the mean free path / of the electrons 
is much smaller than the sample size L, may be described by random matrix theory (RMT). 
The major feature of the spectra, known long ago from RMT is a level repulsion. It 
means that the probability, -P(e), to find two consecutive levels a distance e apart tends 
to zero as e decreases: -P(e) 0. Despite great progress in the theory, closed analytical 
expressions for the distribution function (DF), P(e), are not available (except small e be- 
havior). It was demonstrated, however, in a vast amount of numerical experiments [Q, that 
a simple Wigner surmise, obtained for a 2 x 2 random Hamiltonian, is an excellent approx- 
imation of the -P(e) of large N x N matrices. Originally the Wigner DF was obtained for 
three distinct symmetry classes of the Hamiltonian: Gaussian orthogonal, unitary and sym- 
plectic ensembles ((GOE),(GUE) and (GSE), respcectively) [||. Subsequently the crossover 
ensembles from one pure symmetry type to another were introduced and investigated In 
the present study we are especially interested in a crossover from GOE to GUE, which cor- 
responds to a gradual breaking of the time-reversal symmetry. To this end it was proposed 
to study a random Hamiltonian of the following structure: 

= i/f) + Mf), < a < 1, (1) 

where Hg^^ and are, respectively, symmetric and antisymmetric real random N x N 
matrices. For a = one has GOE, for a = 1 the GUE-case. It was shown recently 
that the time-reversal symmetry breaking parameter a may be uniquely connected to an 
Aharonov-Bohm flux, $ = carried by a flux line, which penetrates the system. The 
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latter, by gauge invariance, may be related to a change of boundary conditions, imposed 
on wave functions in the angular direction (around the flux line): \E'(27r) = \l'(0)e*'^ [0. 
According to Ref. the relation between a and has the form 

where A is the mean level spacing and E^, is the Thouless correlation energy of the system. 

Let us now point out that the first and second derivatives of the energy levels with respect 
to flux (hence with respect to a) have a clear physical interpretation. The first derivative 
(9e„/(9$ is exactly a single level current, carried by each energy level in the absence of time 
reversal symmetry (a > 0). These single level currents manifest themselves, for example, 
in persistent currents through a mesoscopic ring. As was conjectured by Thouless and 
subsequently discussed by Akkermans and Montambaux , a typical single level curvature 
(the second derivative with respect to the flux) may be considered as a measure of the 
correlation energy, hence it is directly proportional to the dissipative conductance of the 
system. The Thouless conjecture is usually written as 

1/2 

(3) 

where the angular brackets (. . .) denote an averaging over a random ensemble. It certainly 
makes sense to ask the question, what are the distributions of the single level current and 
curvature, as a functions of flux, in the framework of RMT. Today one knows only the 
lowest correlators of these quantities like an average value, variance |[T0|JTT[1 or two point pair 
correlator (see e.g. Ref. 10). This knowledge comes mainly from diagrammatic perturbation 
theory, which is usually not applicable to small values of flux. 

The aim of the present article is to derive exact analytical results for the single level cur- 
rent and curvature DF's, using RMT. We shall also find all the moments of these quantities 
and compare them with available perturbative results. Finally we shall calculate the two 
point current-current correlation function for a small flux. As was done in the case of the 
level spacing distribution, we shall use a simple but exactly solvable 2x2 random matrix 

3 




model. We prove, however, that some key features of our results are not restricted to this 
toy model but apply to N x N random matrices as well. 

The most surprising features resulting from the present study are the following: Let us 
define a single level curvature to be proportional to the second flux derivative of energy 



c ~ 



(numerical factors will follow in the body of the article). Then the flux dependent distribu- 
tion function of the curvature has the form (for <^ 0^ = y A/ (S-^c) ^ 1) 



-1/2. 



c< 1, 



(4) 



This distribution implies the following flux dependence of the moments of the curvature 



Ec', 

E-(0,/0)2-™; 



m = 1, 
m = 2, 
m > 3, 



(5) 



As is shown in Appendix A, this structure of the moments is not an artefact of the simple 
model, but may be rigorously derived from a general N x N random Hamiltonian, Eq. (^. 
Based on this, one may suggest the universality of the large curvature, |(9^e„,/(90^| ^ Ec, 
behavior of the distribution, Eq. (^). In contrast, the small curvature part, |9^e„/c?0^| <C Ec, 
may be a specific property of our toy model. 

Equation (||) definitely contradicts the frequently used expression, Eq. (^, when eval- 
uating the average over the disorder as arithmetical mean. Indeed, according to Eq. (^ 
the second and higher moments of the curvature diverge at zero flux (or, in other words, 
in the GOE ensemble). This is a consequence of the absence of the exponential tail in the 
distribution function at exactly zero flux. In this case the behavior Po{c) ~ c~'^ continues 
up to infinity. As was already mentioned, such a behavior is a common feature of GOE 
and not a result of an oversimplified model (see Appendix A). Thus when evaluating the 



Thouless energy, Eq. (|), one should either use a geometrical mean, as already pointed out 
by Thouless [^, who assumed a simple Lorentzian distribution of the curvatures, or use 
another measure of the sensivity of the spectrum to variations of the boundary conditions: 









4>=o 



The same divergence affects dramatically the current-current correlation function. 
Namely, it will be shown that for 0, 0' -C (pc one has 

Differentiating this result with respect to and 0' and then putting 0' = 0, one returns back 
to the just quoted second moment of the curvature. This expression is absolutely unexpected 
from the point of view of perturbation theory. The latter assumes rather E^(f)(j)' (without 
logarithm) for the above defined correlation function. It would be definitely interesting to 
see if such a behavior exists in the framework of supersymmetric calculations. In fact, these 
calculations have already been done in Ref. ||I6[, but, as far as we know, only in GUE, where 
we do not expect anything unusual. 

Curvature distributions have already been investigated in a number of works |]T3|,|T^,|T3[ , 
however, in a very different context. The considered Hamiltonians had the structure H{X) = 
Hi + XH2, where H{X) belongs to the same universality class in the whole range of A ||15|| . 
The curvature is defined as the second derivative of energy with respect to the parameter A. 
Subsequently, curvature distributions for three pure symmetry classes were studied. Since 
the above defined parameter A does not break the time-reversal symmetry, it can not play 
the role of a magnetic flux. Thus, it has no meaning, for example, to look for a current 
(defined as a derivative with respect to A) distribution. For reasons, which become clear in 



Appendix A, our results for a = coincide exactly with the GOE results of Refs. |]I3|,|T^,|T5[ . 
We recover the tail of the curvature distribution, first discovered in Ref. WM, and Eq. 



of the present work may be found in Ref. However, for any a 7^ our conclusions 



are very different from those of Refs. |l3|JT^Jl5|] . For example, in GUE {a = 1) we have find 



a Gaussian tail of the curvature (with respect to flux) distribution. At the same time the 
distribution function of the curvatures, defined with respect to the parameter A, decays only 



as the fourth power (see Appendix A and Ref. p^). 

The present article has the following structure. In section y we specify the 2x2 model 
based on Eq. (|I]) and re-derive the known results concerning the energy spacing distribution. 
In sections |TTl| and |V| single level current and curvature DF's are derived. Finally in section 
we discuss the possible implications of the simple model to real physical systems. In 
Appendix A the moments of the curvature for a general N x N random matrix model are 
considered. A summary of the results for pure ensembles (GOE and GUE) is given in 
Appendix B. 



II. THE MODEL AND ENERGY SPACING DISTRIBUTION 



Consider a model based on H^'^\a), as defined by Eq. ([T|) 



a) 



V 



Xi + X2 X3 + laxi 

X3 — iaxi Xi — X2 

where Xj (j = 1 . . . 4) are real random variables, with a Gaussian distribution law 

1 



(6) 



P{x,) 



X, 



exp 



(7) 



The variance of the distribution f ^ and the time-reversal symmetry breaking parameter a 
are the two free parameters of the model. In the end they should be related to physical 
observables such as mean level spacing and magnetic flux. Let us, however, postpone this 
discussion until section 0. The spectrum of the Hamiltonian, Eq. (^, is given by 



1/2 



e± = xi±. {x\ + 0:3 + a^x^"^ ' , (8) 

1 /2 

and the energy spacing e(a, x) by e(a, x) = e+ — e„ = 2 {x\ + X3 + 0:^X4) . Let us consider 



an energy spacing DP 



r 4 

-Pa(e) = / (^(e ~ e(ct)X)) P{xj)dxj. 



(9) 



The calculation of the integrals is straightforward, finally one obtains 
PJe] 



exp 



erf 



In two limiting cases one returns again to the familiar distributions: for a = 



(10) 



Pnfel 



4y2 



exp 



(GOE), 



(11) 



and for a = 1 



4v27rt;3 



exp 



(GUE), 



(12) 



It is a well-known fact that in the case of GOE, Po{e) oc e for e <^ v, whereas in the absence 
of time-reversal symmetry (GUE) the level repulsion is stronger: Pi(e) oc e^; {e <^ v). In 
the intermediate region < a < 1 one has 

va 

(13) 



1 / 

Pa{c) ^ - exp - — 



e y 1 _ 



e > 



va 



I \2vJ 2^1^ 

The level repulsion is quadratic for small energy intervals, and becomes linear for larger 
ones. The moments of the distribution are given by 



(-1)"2 



nr)37i+ 



1 9^ 



t ""^/^ arctan - 



'1 -a^ 



|t=i ; m = 2ra - 1, 
; m = 2n. 



(14) 



In particular for the first moment one has [17 



a + [1 — a^) ""^/^ arctan ' 



(15) 



This is a smooth monotonous function of a, which varies from (e)o = y f f up to (e)i = 
v^J^2. The fact that it is almost constant will be useful in section 0, where we shall try to 
give a physical interpretation of the results. 



7 



III. SINGLE LEVEL CURRENT DISTRIBUTION FUNCTION 



Define now the single level currents in a 2 x 2 model as 

«±(a,x) = — = ±- r. (16) 

oa e(a;,xj 

We will look for a DF, Pa{i) of i = i+ > 0. Obviously the distribution of i- is the same 
(up to the minus sign of the argument). If one is interested in the DF of z, which may have 
either sign, one simply has Pa{i) = -Pq(K|)/2 (the coefficient 1/2 takes care of the correct 
normalization). To evaluate Pa{i) let us first calculate the joint (current and energy space) 
distribution, Paii^e). Besides technical advantages, this way of calculation provides some 
additional information. Namely, one will be able to identify those energy spacings e, that 
are mainly responsible for a given current, i. The joint DF is defined as 

/oo ^ 
5{i — i{a, x))(5(e — e(a, x)) P(xj)dxj. (17) 

3=1 

After some calculations one gets 



2 (t; 72^)3 y/2m V 



P^{i, e) = Q{e - 2m)-j-^=^^= exp ( . (18) 



Here d{x) is a usual step function (remember that both i,e > 0). To get the current 
distribution, Pa{i), one should, according to Eq. ([T^), integrate the last expression over e. 

P„(0 = i^=^(2m)2 r(l + t)3/^ exp (~^itW + til + a^) + 1)) dt, (19) 
2 (t>v27r)'^ Jo \ 2v^ J 

where a variable t was introduced as e = 2ai{t+ 1). The last integral is not known in special 
functions, except in the two limiting cases a = and a = 1 (see below). However one can 
work out it's asymptotic behavior in various regions 

r r(i/4) 



29/4 Via' 



i -C va; (e ~ v) 



P^{i) ^ - <{ ( -j ; < i < (e ^ m) (20) 

^^exp[-^); {e^av) 
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The values of energy spacings e, which provide the main contribution in each case are des- 



ignated in brackets. Equation (|20| ) shows that the single level current DP has an integrable 
square root singularity at small currents. Realizations with energy spacings of the order 
of an average one are mostly responsible for this singularity. In the intermediate region, 
which exists only if a <^ 1, the current DP decreases as Pinally, for large currents the 
distribution has a Gaussian tail. Por a <^ 1 the last two regions arise due to realizations 
with extremely small energy spacings. 

Por the two pure cases one can calculate the distributions analytically: for a = 

PoW = -26 (-) , (GOE) (21) 



V \v, 

which is evident without any calculations, and for a = 1 



1 1 V2v^ 5 



I 



^'W^^TSVV^^i-iTSj j- («UE, (22) 

where r(a, x) is an incomplete gamma function. The asymptotic behavior of Pi{i) is given 
by the first and the third lines of Eq. (|20|) (with a = 1). The moments of the current 
distribution are given by 



^'n jla ; m = 1, 



-3a In a; m = 2, (23) 

2W2+1 2p/^ + 3\^_ . „ 

-a 1 ; m > 3, 



2™' 

V 2 / m 
for a -C 1 in leading order in a, and 

in , ^ ' (24) 

^"-1-^^ [ 2 ) m + 1/2' ^ 

for a 1 in leading order in (1 — a^). Let us mention particularly that the second moment 

of a single level current is given by (i^)a<i ~ — 3f^Q;^lna. This can be seen directly from 

Eq. (|20|). The leading term comes from the intermediate region of the currents {va <^ 

i <^ v). The main contribution to arises from realizations with very small spacings 

e ~ a^v <^ (e). Starting from the third one, all the moments are determined by a Gaussian 

tail of the distribution {v whereas the average current comes from the opposite 

region i ^v. 



IV. CURVATURE DISTRIBUTION FUNCTION 



Following the same scheme as in the previous section, consider now the single level 
curvature DF. 



± 



e a, X 



e^(a, x) 



(25) 



Again we consider a distribution of a positive defined quantity c = c+ > (it is indeed 
positive as e > 2ax4, according to Eq. (||)). If one wants to include also c_, one should 
again make the distribution symmetrical with a proper normalization. The joint (curvature, 
energy spacing) DF, Pa{c, e), may be found in elementary functions in the closed form 



Pa{c, e) 



TC 



,3/2 



e — sa c 



4 {vV^tt)'^ ^/c 



exp 



VTTb 



exp 



1 



a 



8f2 2a2 



B 



exp 



1 



a 



8v^ 2a2 



B 



(26) 



where B = yl — 8a^c/e. One now integrates the last expression over e and obtains the 
following asymptotics for the curvature DF: 



f r(i/4) ru 

29/4 V c 



Pa(c) 



1 

— < 

V 



3(- 



25/2^4 C 



- exp 




c <^ f , 



f < c < -, 
a 

f V 

- <c< — ^ 

f 

a 1 — 



v) 
av) 



(27) 



where Qa = (Ta^ + ISa^ + 7)/(l + a^). The first two lines in Eq. ( pT]) look very similar to 
those of the current distribution (cf. Eq. (|20|)). The reason is that for c <^ f/a the last 
term in the expression for the curvature, Eq. (pS]), may be omitted, thus one has a trivial 
relationship between current and curvature, i = ac. The large curvature {c ^ v/a) tail of 
the distribution is affected by the last term of Eq. (|25|) , resulting in a complicated form of 
the tail of the DF, Eq. (|27|) . For a = (GOE) the tail disappears completely, in this case 
one can get an exact result for the curvature distribution (see also Ref. [jl5|) 



3 2f fl / c 



D-5/2 f — 



2v 



(GOE) 



(28) 
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where Da{x) is a Whittaker parabolic cylinder function. The asymptotic behavior of this 
DF is given by the first two lines of the general expression, Eq. ([27|). Due to the absence of 
an exponential tail in a GOE curvature distribution, Eq. (^), all the moments of it, starting 
from the second one, diverge. Indeed, a direct evaluation of the moments in leading order 
in a results in 



—3 In a; 

2m/2+l 



m + 3 



E ^ 



m = 1, 
m = 2, 

m > 3, 



(29) 



for a < 1 (cf. with Eq. (^), and 



2^/2 



m + 3\ g 



p=0 



m + 1/2 + 



(30) 



for a 1 (cf. with Eq. (|2^) ); I p j is a binomial coefficient. As in the case of the moments 
of the current, all the moments of the curvature, starting from the second one arise from 
realizations with small energy spacings e ^ (e) (if a <^ 1). At a = only the first moment 
exists, whereas all higher moments diverge. This unexpected fact will be discussed in more 
detail in the next section. 



V. DISCUSSION OF THE RESULTS 

As already mentioned in section |, the Wigner surmise obtained for 2x2 matrices, 
works extremely well also for a large A^. To establish this connection one should relate 
the phenomenological parameter f ^ - the variance of the distribution - to an average level 
spacing A. One simply demands that 

{e)a = A. 

Strictly speaking the average spacing (e)^ is a function of a (magnetic fiux), although 
the mean level spacing A is presumably a constant, independent of external parameters. 
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However, as we noticed after Eq. (|T5D, the dependence on a is very week (especially for 
small a). Using this fact we shall disregard its a-dependence and just admit 



where we have used Eq. (0). Being honest, one should re-identify the parameters for each 
value of a separately. This procedure, although trivial, is not transparent enough for our 
illustrative purposes. 

Having an energy spacing DP as an example, one may hope that the single level current 
and curvature DF's, derived for a 2 x 2 system, may be suitable for larger systems as well. 
To support the last statement let us put forward the following arguments. As we have seen 
in previous sections, all the moments of a single level current and curvature DF's, starting 
from the second one, arise mainly from realizations of a random Hamiltonian with very 
small gaps. For these realizations the 2x2 ansatz is supposed to be essentially correct, 
because for the close pair of levels only their mutual interaction appears to be important. 
In Appendix A we prove that the small flux behavior of the moments is indeed observed in 
the general N x N model as well. The first moment, however, is determined by realizations 
with an energy gap of the order of the average one. In this case the 2x2 scheme need not 
be precise. Thus one should not trust the value of the first moment, but rather connect it 
phenomenologically with the microscopic characteristics of a system. Following Thouless, 
one may relate a typical second derivative (not r.m.s !) with respect to flux at zero flux to 
a correlation energy 







lO 






</)=0 



cf. with Eq. (H). On the other hand, we had (see Eq's. (^), (^)) 




(31) 




Using the definition of a mean level spacing, Eq. (^Tf), one obtains 




(33) 
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This should be compared with the conjecture of Dupuis and Montambaux Eq. (H) (with 
= 2). As one sees, the agreement is extremely good, the slight discrepancy may be 



2x2 model are applicable for larger systems as well. Let us discuss the further consequences 
of this rather strong assumption. 

First of all one notices that the a = 1 (GUE) case corresponds to the value of a flux 



0c = Y^/(2-E'c)- This value is well-known as a correlation flux. Up to this flux a typical 
level may change parabolically, without crossing other levels. At = 0c the first avoiding 
crossing event usually happens, and the simple 2x2 scheme obviously breaks down. It was 
demonstrated numerically 0, that at ~ 0c the crossover to GUE is indeed practically 
completed. This shows that the applicability of a 2 x 2 model for < < 0c (0 < a < 1) 
is quite reasonable, as well as the identification of the = 0c point with GUE. 

Consider now the second moment of the single level current in GUE (or, the same, at 
= 0c). Using Eq. (P^D, one obtains ((i^)Q,=i)^^^ = vJs/b, or in physical parameters (using 



This result is also well-known from perturbation theory (up to the numerical coefficient) 



Being thus convinced that the obtained results lead to reasonable predictions for real 
physical systems, let us discuss the most surprising feature of the considered DF's: At zero 
flux all the moments of the curvature, starting from the second one, diverge. Thus, when 
calculating the correlation energy from the curvatures, one has to use another measure for 
their typical value than just the root mean square, like Eq. (p2D or the geometrical mean 
proposed by Thouless 0. 

The consequences of the discussed divergence are, however, deeper than just the necces- 
sity of a more careful definition of the correlation energy. One also should reconsider the 



attributed to the fact that Eq. (Q) was obtained for the large A^ limit. We conjecture thus, 
that with the identifications, Eq's. (|31|), (|33|), the tails of the distributions obtained for a 




Eq's. (0), (0)) 




mm- 
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universal relationship between dissipative and correlation conductances, derived by Akker- 
mans and Montambaux 0. Mathematically this relation was expressed as 

2 



( 



<9e„ 


2 

) = aA( 















) 



1/2 



where bar denotes integration with respect to flux, and a is a universal numerical factor. 
According to the present results this relation can not hold, when the typical curvature is 
calculated as an arithmetical mean. Indeed, using Eq. (|19D and Eq's. (|3l|) , (|33D , one obtains 
for the l.h.s. of the last expression y/ AEc (up to a coefficient of the order of unity), whereas 
the r.h.s. diverges. 

To understand the reason for this phenomena let us consider the two point current- 
current correlation function 



C{a,a') = {i{a)i{a')). 



(34) 



A very similar object was recently considered in Ref. |12|. One can explicitly perform the 
averaging in a 2 x 2 model by integrating over rfx with the corresponding weight, precisely 
as one did in the previous sections. The general answer is cumbersome, but one needs only 
the behavior for small flux. In this case one easily gets 



C(a, a') ~ —3v aa' ln(a + a') ; 



a, a' -C 1. 



(35) 



Putting here a = a', one returns again to the expression for the second moment of the 
current (the second line in Eq. (|23D ). On the other hand, differentiating Eq. ( ^Sf ) with 
respect to a and a' and then putting a = a', one recognizes the second line of Eq. (^9]). In 
physical parameters Eq. (|35|) may be rewritten as 



< 0c- 



(36) 



This should be compared with the corresponding result of the perturbative calculations 

2 



-2 772 - -/ 



A 

,7, 



<0c, 



(37) 
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where 7 is a cut off in perturbation theory, which is usually supposed to be of the order of 
A 0. The discrepancy between the two results is rather dramatic. Whereas Eq. (^) leads 
to a finite second moment of the curvature (~ E"^), Eq. ( p6D results in a divergent second 
moment. Let us also point out that the perturbative result obviously may be expressed 
in a form C(0, 0') = /(0 + 0') — /(0 — 0'), which may be traced back to Diffuson and 
Cooperon channels in the diagrammatic expansion. The Eq. (|36D does not allow such a 
decomposition. This might be a point where the present scheme contradicts the derivation 
of Ref. [^. Indeed, it was assumed explicitly that the Diffuson-Cooperon decomposition 
(which is certainly correct for a large flux ^ 0c) is also valid in the vicinity of zero flux. 
According to the present consideration this is not the case. 

It is not clear at the moment whether the discussed divergence has a real physical mean- 
ing, but if so, it might cause difficulties in numerical calculations of correlation functions in 
GOE. We conclude that further analytical (both RMT-like and supersymmetric) calcula- 
tions and numerical work are necessary to clarify this unexpectedly controversial issue. 
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APPENDIX A: MOMENTS OF THE CURVATURE IN A iV x MODEL 

Let us show now that the flux dependence of the moments of curvature, derived for a 
2x2 model, Eq. (|29|), and summarized in Eq. (|^), may be obtained from a general N x N 
model. Consider an x random matrix Hamiltonian, given by Eq. (0). Without loss of 
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generality one may assume that its spectrum is not degenerate. Then, using second order 
perturbation theory, one obtains the following exact relationship 



^^=2eJ — ; , , (Ai) 

where ekicn) and \k,a > are eigenvalues and eigenfunctions of the full Hamiltonian, H^^\a). 
As it is well-known from RMT statistics of eigenvalues and statistics of eigenfunctions 
are completely independent of each other. 

Let us first consider the case of exactly zero flux {a = 0, GOE). In this case the energies 



in the denominator on the r.h.s of Eq. ( |Al| ) are eigenvalues of Hg^\ whereas in the numer 



ator one has matrix elements of H^^\ Hence, matrix elements. 



< 0,k\Hf^\n,0> 



and 



eigenvalues, = efc(O), may be considered as independent random variables. The statistic 
of the eigenvalues is given by a Wigner-Dyson (GOE) distribution |Q 

PAr(ei,...,eAr) =const X exp (-^^e^J [| |e„-efc|, (A2) 

\ k=l ) l<k<n<N 

whereas an exact form of the matrix element distribution is not important for our pur- 
poses. One is now in a position to consider the moments of the random variable 
c = |(9^e„(a)/(9a^|a=o|- Doing this, one will be interested only in the manner of divergence, 
omitting all the prefactors as well as less divergent terms. Raising Eq. ( |Al| ) to the m*^ power 
and averaging, one obtains 

in-N^' [[^^^de,de, + ..., (A3) 
J J ei - e2 p 



where averaging over matrix elements. 



< 0,fc|M^V,0 > 



leads to some omitted constant 



prefactor, and ". . ." denotes less divergent terms, arising from the non-diagonal contribu- 



tions, like 



-R3(ei, £2, es) J J J 
aeiae2de2 



1^1 ~ ^2 1 1^2 — ^3 



etc. Here Ri{ei, . . . , e^) is an i-point correlation function [Q, for example 

roc /"OO 
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In the limit of large A^, the correlation functions Ri depend only on differences of the 
eigenvalues. Then the integral over ei + €2 in Eq. (^^) leads to some constant of the order 
of A^, and one finally obtains 

- /f ^de, (A5) 

where e = |ei — 62!. Taking into account the well-known GOE result, -R2(e) ~ e (e ^ A), 
one notices that all the moments of the curvature, starting from the second one, diverge. 
This is in exact agreement with the result for a 2 x 2 model, but now one has demonstrated 
the validity of this statement for a general N x N model. 

The above developed scheme is applicable without any changes for the case where a 
is not a time-reversal symmetry breaking parameter, hence H^'^\a) belongs to the same 
universality class in the whole range of a [|T^,|T3|. Indeed, in this case the matrix elements 



of the perturbation (numerator on r.h.s of Eq. and the eigenenergies (denominator) 



can always be considered as independent random variables. One can, for example, repeat 
exactly the same arguments for GUE with a single modification: R2{^) ~ e^- Then one 
obtains the divergence of all the moments of the curvature starting from the third one. This 
perfectly agrees with the result for the tail of the curvature distribution, proved in Refs. 
|1^,[15|: P(c) ~ c""^. In the same way for GSE, using -R2(e) ~ e"^, one obtains the divergence 
of the moments starting from the fifth one, in agreement with P{c) ~ |jl3|,|l5|. This is 
the reason, why the result for the GOE curvature distribution, Eq. (pH]), coincides exactly 
with that of Ref. p| . 



However, when a plays the role of a magnetic flux, it does change the symmetry class of 
the Hamiltonian. Then the energies efc(a) are eigenvalues of the sum i/f^+zai/f), and 
thus no longer statistically independent of the matrix elements, < a, k\H^^\n, a >. We are 
going to show that in this case, for a 7^ 0, all the moments converge. This in turn means 
that, contrary to the case of Refs. [Il3| , |l5i , the curvature DP has an exponential tail. For 
small a and small |ei — 62] one can show that 

|ei(«) - e2(«)| ^ ((ei - e^f + a^X^^'^^ 
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where the quantity X depends only on the matrix elements (and not on the energies), hence 
it is statistically independent of the zero-flux energies, e^. Repeating again the steps, leading 
to Eq. ([A5|) one obtains 



const; 171 = 1, 

Ina-i; m = 2, (A6) 
rn> 3, 







where f{X) denotes the (unspecifled) distribution of matrix elements. In Eq. ( |A6|) we used 
again -R2(e) ~ e- Thus we have shown that the small flux behavior of the moments of the 
curvature, initially derived for a 2 x 2 model, Eq. (|), is valid for large N as well. 

APPENDIX B: LIST OF DISTRIBUTIONS FOR PURE ENSEMBLES 



GOE; a = Q 
Spacing distribution: 



Current distribution 



W-^exp , (Bl) 



(e'")o = t;-23™/2r (^^^) . (B2) 



Po(0 = -24-)' (B3) 



V \v 



{no = 0. (B4) 



Curvature distribution 



(c)o = t^y|; (Oo = oo; m>2. (B6) 
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GUE] a = 1 
Spacing distribution: 



e2 / 



23m/2+l /™ I Q 



TT 



Current distribution 



m = - 



1 1 ,/v^j./5 / z 



{i )i = V — 



Vtt V 2 /m+1/2 



Curvature distribution 



( r(i/4) /V 



29/40F V c' 



0r V 2 /^m + l/2+p' 
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